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Abstract 

We use action-angle variables to describe the basic physics of coherent scalar field oscillations 
in the expanding universe. These analytical mechanics methods have some advantages, like the 
identification of adiabatic invariants. As an application, we show some instances of potentials 
leading to equations of state with p < —p/3, thus contributing to the dark energy that causes the 
observed acceleration of the universe. 
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I. INTRODUCTION 



The standard model of cosmology assumes the Friedmann-Robertson- Walker metric 

dr 2 



ds 2 = -dt 2 + R 2 (t) 



+ r 2 dn 2 



hr' A 



(1) 



corresponding to an homogeneous and isotropic universe. Here k is the curvature signature 
and R is the expansion factor, whose time change is given by the Friedmann equation 

We have written the equation in such a way that the cosmological constant is included in 
the total energy density p?- 

There are different contributions to pr- Matter and radiation are among them. They can 
be introduced as a fluid with pressure proportional to the energy density, p — wp; w — 
corresponds to matter and w = 1/3 to radiation. Recent results coming from high-redshift 
supernovae cosmic background radiation Q], and galaxy survey [3] suggest contributions 
with w ~ — 1. This is the so-called dark energy that leads to the present acceleration of 
the universe. This acceleration may be produced by a cosmological constant, since it has 
an equation of state with w = —1. However, there might be other causes, like quintessence 
^], cardassian expansion etc. Of course, there might be several components building up 
dark energy. 

In this paper, we would like to discuss some aspects of the contribution to the energy 
density of coherent scalar field oscillations in a potential. The physics of such oscillations 
was analyzed in a pioneer work by Turner p]. The purpose of our paper is the following. 
First, in SectlHl we find the analytical form of the field oscillations in a couple of interesting 
cases. In Sect JIIIl we treat the general case and derive relevant physical results using the 
action and angle variables of analytical mechanics Q|. We find this language very useful, 
particularly when an adiabatic change of the potential is present. Another objective of 
our paper is to show that scalar field oscillations may contribute to the dark energy of the 
universe. We do this in Sect II V[ where we present some instances of potentials leading to 
an accelerated universe. 



2 



II. ANALYTICAL SOLUTIONS FOR SCALAR FIELD OSCILLATIONS 



A classical scalar field in a potential has the Lagrangian density 

C = -\d,<pd»<p-V(<P). (3) 

In the metric (JTJ), and considering spatially homogeneous configurations, the equation of 
motion for the field is 

4> + 3#0 + V'(4>) = 0. (4) 

Here dot means time- derivative and prime means d/dcf). 

In this section we concentrate on two particular potentials, V ~ <p 2 and V ~ 4 , where 
one can find analytical solutions. While the solution for the harmonic potential is in the 
literature |8|], the solution we find for <ft A is new, as far as we know. 

For the quadratic potential 

V = \m 2 <? (5) 

we have harmonic oscillations, cj){t) = Asm(mt + <p), in the case H = 0. 

When H ^ 0, but in the case H <ti m and H/H -C m, we expect a time dependent 
amplitude and, perhaps, a different frequency. We make the following ansatz 

4>(t) = A{t) sin(Amt + <p). (6) 

We find A = 1 and the following equation 

A = -\ha. (7) 

The energy density p corresponding to evolves as 

p oc A 2 oc R- 3 (8) 

where in the last proportionality we have used the solution of (JTj). This corresponds to the 
behavior of non-relativistic matter in the expanding universe. 
The second potential to be solved is 

V = -m 2 <f)\ (9) 

2 

We shall work in complete analogy with the previous example. We first find the solution to 
(gj) for H = 0, which is 

<f>(t) = B sn(Bmt + cp) (10) 



where sn is the Jacobi elliptic function sn(Bmt + ip, — 1) 
For H we choose an ansatz of the form 

4>(t) = B(t)sn(B(t)Xmt + <p). (11) 

The calculus here is more complicated than in the previous example, when introducing (jll|) 
in (£Q). We obtain a differential equation for the evolution of the amplitude B: 



dt 2 V tj dt 
This equation can be easily solved for the usual form H = ct~ l , where the value of the 
constant c depends on which kind of energy dominates the universe (c = 1/2 for a radiation- 
dominated universe and c = 2/3 for a matter-dominated one). From ()12j) we obtain 



d(Bt) 
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oc t 2 - 3c (13) 



dt 

that gives 

B oc r c . (14) 
Since we are assuming H = ct~ l , which gives R oc t c , we conclude that, for any value of c, 

B(xR-\ (15) 

Also, we get an equation for A, which now depends on the relative variation of B 

A =( 1 + l ( ) =jh < 16) 

where in the last equality we took into account (|14)1. 
The energy density in this case is given by 

p = \Q>y+ l -m\4>Y= l -m*B* (17) 

where we have used properties of the Jacobi elliptic functions [9]. Thus, taking into account 
(JT5J), we have 

p oc B 4 oc R-\ (18) 

This type of i?-dependence corresponds to a fluid of relativistic particles, or equivalently, to 
radiation. The form of the oscillations in a 4 potential with small friction H was also found 
in Q], in the context of preheating after inflation. The method used in was to make 



a conformal transformation of the space-time metric and the fields, as well as a number 
of re-scalings. Even if we get the same final results, we present our method because it is 
somewhat simpler. 

The two examples discussed here are particular cases in which it is possible to analytically 
solve the equation of motion that describes the oscillations of the field in the expanding 
universe. For other potentials that can be considered, it might be impossible to find an- 
alytical solutions to (@J), so the method applied before might not succeed in obtaining the 
dependence of the energy density p on the scale factor R. For this purpose, one needs to 
find another method, in which the evolution of p(R) can be obtained without necessity of 
solving the equation of motion (jJJ). This is what we develop in the next section. 



III. ACTION- ANGLE FORMALISM 



We are concerned with the oscillations of the 0-field about some minimum of the potential, 
but we shall not restrict to the case that the oscillation amplitude is small. Thus, we are 
faced with a system with a periodic motion, whose details might be complicated. Often, we 
are not specially interested in these details, but rather in the frequencies of the oscillations, 
the contribution to the energy density, etc. A method of analytical mechanics tailored for 
such a situation is provided by the use of action-angle variables. 

As before, it is convenient to start with H = in (J3J); we have then a conservative system 
with hamiltonian density H. 

ft(0,n) = in 2 + v^) = P (19) 

with the momentum II = <fi and energy density p. When having a periodic motion, one 
introduces [?J the action variable 

J = ^Ud(p (20) 

where the integration is over a complete period of oscillation. Using II = \/2(p — V) from 
(fT9*|) . J can be written as 

J = 2 v / 2(p - V)d(j> (21) 

where mi „ and (p ma x are the return points, V(0 min ) = V((f> max ) = p. J is chosen as the new 
(conserved) momentum in the integration of the Hamilton- Jacobi equation. The generating 



5 



function given by the abbreviated action 

W = I Ud(p (22) 



allows to canonically transform (0, II) into (a, J), with the angle variable defined by 

a = jj. (23) 

Since W does not explicitly depend on time, the new hamiltonian 7i coincides with the old 
one Ti — p. The energy p is only a function of J, which amounts to say that a is cyclic and 
J is constant 

j-f a -, (24) 



The other Hamilton equation is 



■ 9H(J) 

«=-£ 1 = ->, P5) 



with v — i/( J) constant. We can integrate to obtain 

a{t) = vt + a(0) (26) 

and in a complete period, a(r) = a(0). In this way, we identify v in ()25j) with the frequency 
of the motion: 

Until here, we have reminded ourselves of the standard approach that uses action-angle 
variables to find the frequency [3]. 

The realistic case of the expanding universe, with H ^ 0, can be treated with the same 
method, provided H and H/H are small compared to the frequency of the oscillations, 
that is to say, 

H^u- HjH^v. (28) 

Notice that for the usual form H = ct _1 , H/H ~ H so that the two conditions in (J28|) are 
actually the same. The conditions (J28j) ensure that the motion is almost periodic, and that 
it makes sense averaging over one cycle. The if-term in (0J) represents a time- dependent 
friction, so we expect energy to decrease; however, we are assuming this friction small enough 
so that we can consider the energy constant in one-cycle period. 
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We still define J as in (j2~Uj) and (}2"T]) . and also the energy density p as in (fT§|) . The 

decrease rate of p can be obtained rewriting the equation of motion 

t - " 3 ^ 2 < 29 > 

When averaging the r.h.s. over one cycle, we can take H out of the average, (H<p 2 ) ~ H((p 2 ). 
The average of <fi 2 is related to J: 

(this is essentially the virial theorem). Thus, 

(?) = p. (3D 

With the help of reads 

What is remarkable about this equation is that the change in the action variable J is 
independent of the form and details of V. The equation can be readily integrated to show 
that J dilutes volume in an expanding universe 

J oc R' 3 . (33) 

Eqs. fl32)l and ()33|) are one of our main results. To find more, we notice that the pressure can 
also be averaged in one oscillation period, 

P = <^ 2 > - (V) = (<j> 2 ) -p = wp. (34) 
To calculate w, we first use (|3T|) and get 



W = -p-dJJdp~ l (35) 



and, introducing (|21|) . we finally obtain 



1. (36) 



ffZZip-v)- 1 **!' 

Eq. (|36|) coincides with the corresponding expression in Ref . || , where another formalism was 
used. 
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There is another useful relation that can be easily obtained using our formalism. When- 
ever w is constant, we can integrate to obtain 

p oc J w+l . (37) 

Then, making use of (}3*3*j) . we get the well-known relation 

pocR- 3{w+1) . (38) 

As an example, let us apply the method described above for the following power-law 
potential 

V = a<f) n (39) 
for which the action variable can be exactly calculated 



J 



n 



+ 2)r(| + i) 

Using (|2*7j) . we get the frequency of the motion 



2V2^r(i)" v« 
and using (J35j) . we get the parameter u> 

n — 2 , , s 

w = . 42 

n + 2 K ' 

Let us show that these results coincide with what we obtained in Sect|H] For n = 2, we get 
from (pf2*|) that w = 0, which corresponds to the equation of state for non-relativistic matter, 
p = 0. Equivalently, from eq.flHBJ), we obtain that 

p oc iT 3 (43) 

which is the same as (JBJ). Still for the harmonic potential, n = 2, taking a = |m 2 , we can 
calculate J by making use of (f47Ij) 

J = 2tt— (44) 
m 

which is nothing else than the number density of particles of mass m, except the factor of 2tt. 
We now see that, in this case, (J3l?j) is equivalent to assert that the number density decreases 
as R~ 3 in an expanding universe. 
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For n = 4, (}4*2j) gives w = 1/3, and (J58|) gives p oc i?~ 4 which coincides with (fTHj) and 
corresponds to the equation of state for relativistic matter. 

As expected, we have recovered the same results of Sect |TT] without the need to solve the 
equation of motion, eq(j3J). For this reason, the action-angle formalism will allow us to study 
more general cases of potentials for which there is no analytical solution to (@J) . This is what 
we will do in the next section. 

As another application of the action-angle variables, we discuss the issue of adiabatic 
invariants. When a parameter a of the potential changes with time slowly compared with 
the natural frequency of the system, 

- < v, (45) 
a 

we say it changes adiabatically. When there is such an explicit time change in the potential, 
W in depends on time, and "H and p do not coincide any longer 

dW 

H = p + — = p + W, a d (46) 
where the derivative W A = dW/da is taken at constant J. Instead of we now have 

j = -3HJ - (dn/da) 

V 1 1 (47) 

= -3HJ- ({dWjda)a). 

When averaging this identity, due to (j43|) . we can consider d as constant during a cycle, i.e., 
((dW t a/da)a) ~ {{dW A / da))a. The average (dW ta /da) vanishes, since W yCL is a periodic 
function of a. It follows that 

j = -3HJ, (48) 

namely, eqs. ()32|) and ()33|) are still valid, although we stress that V is changing with time. 
For H = 0, the system is conservative and we have J = 0, i.e., J is identified with an 
adiabatic invariant Q|. 

An instance where we can apply these results is the invisible axion model. In this model 
we have a time-dependent axion mass, m = m a (t) in (jSJ). Even if p(t) and m a (t) are 
complicated functions of time, we conclude after our study that, provided the change is 
adiabatic, the combination 

oc R- 3 (49) 



2tt m a (t) 

behaves as the axion number density. This behavior is what is expected on general grounds. 
This result was derived in Q when calculating the relic axion density. What we have done 



in the present paper is a kind of generalization of (|49jh that may be applied when having 
other potentials and/or other types of adiabatic change. 



IV. DARK ENERGY FROM FIELD OSCILLATIONS 

In the standard model of cosmology, besides eq.(J2J), there is another equation that one 
can get from Einstein equations. It can be put in the form 

R AttG , . , , 

^ = - — (P + 3p). (50) 

For the simple equation of state p = wp, we see that w < —1/3 leads to an accelerated 
universe. We will show in this section that such values of w can be obtained from an 
oscillating scalar field whose energy density dominates. An application based on this idea 
in the context of inflation was developed by Damour and Mukhanov in (see also jl^|). 

In order to check whether a given potential may lead to acceleration, we will calculate 
w using the expression (|35p. that involves an average over one cicle |l6j . The value of w 
depends on the form of the potential and on the return values (f) max and (f) m i n . In turn, these 
change with time due to the friction produced by the expansion of the universe, so that, in 
general, w depends on time. 

A remarkable exception to this dependence on the return values is given by the power- 
law potential (}3*9"j) . As shown in (|42|) . w only depends on the power n. For integer n, we 
do not have oscillations neither for n = 0, nor for n = 1. The potential V = acj) 2 leads 
to oscillations that correspond to non-relativistic matter with w — 0. The next potential 
leading to oscillations is V — a</> 4 , that we saw it implies w — 1/3. Higher values of the 
power n lead to higher values for w. No integer value of n, leading to oscillations, can 
correspond to a fluid with w < —1/3 

In order to have smaller w values, one clearly needs a field that, even oscillating, has 
a slow-roll for enough time. To illustrate that this is possible, we shall investigate a few 
potentials. In all of them, we shall put the center of the oscillations at = 0. In addition, 
we will work with symmetric potentials, so that —<pmin = 4>max = <Po- 

Let us start with the "Mexican hat" potential shown in Figd 

Vi{4>) = Ai(0 2 - vl) 2 (51) 



10 



with Ai a coupling and V\ an energy scale. We consider O > V^Vi, to have oscillations 
around = 0. The parameter w, calculated numerically, is a function of <fio/vi and it is 
displayed in Fig|21 We see that when O is close enough to V2vi, we have values w < —1/3. 
However, an unsatisfactory feature is that one has to tune the return value 0o- 
A potential that does not have this fine-tuning problem is given by 

W> = (52) 

(see Fig|3J). Here, t>2 and v' 2 are energy scales. The values of w are shown in FigHJ where 
we see that w < — 1/3 for 0o > 1.2 f 2. To get the observed acceleration of the universe, we 
should have v ' 2 — 2 x 10~ 3 eV, as expected. Notice that in the past, field oscillations give 
w ~ — 1 always. As happens in the case of a cosmological constant, the contribution we 
are discussing was subdominant in the past, as soon as non-relativistic matter enters in the 
stage and dominates. 

One could find other potentials giving acceleration. One such example was discussed in 



ll(. There, the objective was to get first slow-roll inflation and, after, oscillations provoking 
further inflation. In contrast to [ll|, we work in the period of dark energy domination and 
we would like to know the future of the universe. We see in Fig0]that the oscillations will 
end up giving w = 0, i.e., being harmonic and consequently matter dominated. 

One may think that for all potentials, when the field is close enough to the minimum, 
one ends with the usual harmonic potential (0) that well approximates oscillations around 
the minimum, so that one obtains w = 0. We would like to point out that there is at least 
one exception. Consider the potential 



VM = v> 3 4 e-^ 2 (53) 

with t>3 and v' 3 some energy scales (see FigEJ) - An inspection to the potential shows that 
for O large enough, we should have w < —1/3. We have calculated w and our results are 
shown in FigEl Indeed, we see that w < —1/3 for O > 1-8 v%. For smaller values of 0o 5 
higher values of w are obtained. We notice that in the limit <po/v^ — > the value w = is 
not obtained, but rather w — 1. This corresponds to a fluid with p = p. The reason for this 
very particular behavior is the well-known fact that V^{(f)) in (|53|) can not be developed in a 
Taylor- Maclaurin series around = 0. Again we should check that V3 is consistent with the 
observational constraints. With ^~2x 10 -3 eV the oscillations fit the acceleration. Also, 
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for early times w ~ — 1 so that the oscillations do not contradict the observed evolution of 
the universe. 

Let us finally address the requirement that dark energy should have very homogeneous 
density, with spatial irregularities rearranging at an effective relativistic speed of sound c s , 
given by 

In the main spirit of our paper, we should point out that (|54|) can be calculated following the 
technique of action-angle variables. This allows to find c s without reference to the details of 
the solution to the equation of motion of the field. The way we do it is through the equation 

^2 = d(wp) = pdw/d(j) + w= Vodw/d(p + ^ 
dp dp/d(pQ dVo/dcpQ 

where Vq = V(cj)o) = p and w was found in (}3l>j) or (J2SJ). 

As a check, we have calculated the integrals (|55|) for the potential (0) and obtained c 2 ~ 0, 
as expected. The implications of a scalar field with the potential (0) when it dominates the 
energy density of the universe have been fully investigated in ref. [15[. Here we do not 
pretend to do such a complete study for the potentials of Sect Jl VI since our emphasis is in 
the techniques developed in Sect JIIII Let us point out, however, that in fact we find that 
for the potentials ()51|). (j52j) and ()53|). c 2 is near zero or even negative, which indicates that 
collapse at small scales would happen or even would lead to dangerous instabilities. 

It is not difficult, however, to find potentials giving w < — 1/3 and c 2 ~ 1. One example 

is 

VM = vt-^ 2 + (56) 

We have verified that w from ()36|) and (? s from (jo3j) for the potential V4 in (J51)j) have the 
expected values for a dark energy contribution, for a suitable range for the potential free 
parameters. 
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FIG. 1: "Mexican hat" potential Vi denned in (f5T|). 
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FIG. 2: w as a function of the amplitude 4>o of the oscillations in the potential V\. The dashed line 
indicates w < —1/3. 
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FIG. 3: The shape of the potential V2 defined in (|52|). 
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FIG. 4: w as a function of the amplitude 4>o of the oscillations in the potential Vi- The dashed line 
indicates w < —1/3. 
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FIG. 5: The shape of the potential V3 defined in (|53|) . 
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FIG. 6: w as a function of the amplitude 4>o of the oscillations in the potential V3. The dashed line 
indicates w < —1/3. 
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